Charged Brownian particle in a magnetic field 

Radoslaw Czopnik 
Institute of Theoretical Physics, University of Wroclaw, 
PL-50 205 Wroclaw, Poland 
and 

Piotr Garbaczewski 
Institute of Physics, Pedagogical University, 
PL-65 069 Zielona Gora, Poland 

February 1, 2008 



Abstract 

We derive explicit forms of Markovian transition probability densities for the velocity space 
and phase-space Brownian motion of a charged particle in a constant magnetic field. 

1 Motivation 



An old-fashioned problem of the Brownian motion of a charged particle in a constant magnetic field 
has originated from studies of the diffusion of plasma across a magnetic field jjj], [|J and nowadays, 
together with a free Brownian motion example, stands for a textbook illustration of how transport 
and auto-correlation functions should be computed in generic situations governed by the Langevin 
equation cf. || but also [Q, ||. To our knowledge, except for the paper (|] no attempt was made 
in the literature to give a complete characterization of the pertinent stochastic process. However a 
striking peculiarity of Ref. Q is that the Brownian motion in a magnetic field is there described 
in terms of operator-valued (matrix-valued functions) probability distributions that involve fractional 
powers of matrices. In consequence, we have no clean relationship with the standard formalism of 
Kramers-Smoluchowski equations, nor ways to stay in conformity with the standard wisdom about 
probabilistic procedures valid in case of the free Brownian motion (Ornstein-Uhlenbeck process), cf. 
@! @' Hi- Therefore, we address an issue of the Brownian motion of a charged particle in a magnetic 
field anew, to unravel its features of a fully-fledged stochastic diffusion process. 

2 Velocity-space diffusion process 

The standard analysis of the Brownian motion of a free particle employs the Langevin equation 
— — (3lf + A (t) where it denotes the velocity of the particle and the influence of the surrounding 
medium on the motion (random acceleration) of the particle is modeled by means of two independent 
contributions. A systematic part —j3ut represents a dynamical friction. The remaining fluctuating 
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part A (t) is supposed to display a statistics of the familiar white noise: (i) A (t) is independent of 
it, (ii) (Ai (s)} = and (A* (s) Aj (s 1 )} = 2qSijS (s — s 1 ) for i,j = 1, 2, 3, where g = is a physical 

parameter. The Ornstein-Uhlenbeck stochastic process comes out on that conceptual basis. 

The linear friction model can be adopted to the case of diffusion of charged particles in the presence 
of a constant magnetic field which acts upon particles via the Lorentz force. The Langevin equation 
for that motion reads: 



^- = -f3lt+^lt + ^(t) (1) 
at mc 

where q e denotes an electric charge of the particle of mass m. 

Let us assume for simplicity that the constant magnetic field B is directed along the z-axis of a 
Cartesian reference frame: ~E> = (0, 0, B) and B = const. In this case Eq. (Q) takes the form 



dlt — i* / v 

— = -Alt + A (t) (2) 



where 

(13 -u c \ 
A = \ uo c (3 (3) 
V J 

and lo c = denotes the Larmor frequency. Assuming the Langevin equation to be (at least 
formally) solvable, we can infer a probability density P (it, t\lto), t > conditioned by the the initial 
velocity data choice it = ito at t = 0. Physical circumstances of the problem enforce a demand: (i) 

3 

P(H > ,t|H > o) -* S 3 {it - lt ) ast^O and (ii) P(lt,t\lto) -* (^f^r) * ex P (- "LbT ) as t ~* °°- 
A formal solution of Eq. (S) reads: 



It (t) - e- M lto = I e- A( - t - s) A{s)ds . (4) 
By taking into account that 



cos Lu c t sin u> c t 
e -At = e -0t\ - s i nLUc t cosu c t \=e~ pt U(t) (5) 
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we can rewrite ([|) as follows 



(t)-e- f3t U(t)lto= / e-^ t -^U{t~s)A(s)ds . (6) 



Statistical properties of it (t) — e~ A *T? are identical with those of A (s)ds. In consequence, 
the problem of deducing a probability density P (it ,t\lto) is equivalent to deriving the probability 
distribution of the random vector 
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7?= f i/j(s)1£ (s)ds (7) 
Jo 

where tp (s) = e'^'^ = e'^-^U (t - a). 

The white noise term A (s) in view of the integration with respect to time is amenable to a more 
rigorous analysis that invokes the Wiener process increments and their statistics, |^|. Let us divide 
the time integration interval into a large number of small subintervals At. We adjust them suitably 
to assure that effectively ip (s) is constant on each subinterval (jAt, (j + 1) At) and equal if> (jAt). As 
a result we obtain the expression 



JV-L r(j+l)At 

^=]TV(iAt)/ A {a) da. (8) 

J jAt 



3=0 



Here i? (At) = J^" 1 ^' ~A* (s) ds stands for the above-mentioned Wiener process increment. Phys- 
ically, B (At) represents the net acceleration which a Brownian particle may suffer (in fact accumu- 
lates) during an interval of time At. Equation (0) becomes 



N-l 



N-l 



S^=5>(jAi)^(Ai) = £ 



(9) 



3=0 



3=0 



where we introduce = if) (jAt) i? (At) — xjjjB (At). 

The Wiener process argument 0], allows us to infer the probability distribution of ~s*j. It is 
enough to employ the fact that the distribution of ~S (At) is Gaussian with mean zero and variance 

q 



-B-^B. Then 

m. 1 



i?(At) 



/ I 



\AnqAt 



exp 



B (At) 



4qAt 



(10) 



and in view of ~s*j = ipj B (At) by performing the change of variables in (|T^) we get 



w 



[-?,] = dot [tpj 1 ] w [tfj 



; detVj 



w 



(11) 



det V (s) = e - 3fl(t -^ and ip' 1 (s) = U [- (t - a)] e^- s) we obtain 



1 / \ e M-**)U[- (t-jAt)} ^,\ 2 \ 

■exp — - — — (12) 



AnqAt J e -a/J(t-jAt) ^ I 4g Ai 



and finally 
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Wm ~ [AwqAt e-W-W ) ° XP { AqAte~^-^) ) " (13) 
Clearly, 1?j are mutually independent random variables whose distribution is Gaussian with mean 



zero and variance cr 2 = 2qAte 2/3 '* i At \ Hence, the probability distribution of S = X^Lo* 



again Gaussian with mean zero. Its variance equals the sum of variances of ~s*j- i.e. a = J2j a j 
2qY^i Ate- 213 ^-^^. After taking the limit N -> oo (At -> 0) we arrive at 
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o* = 2q[ dse-W*-') = — (l - e^) . (14) 
Jo m 



Because of 5* = Tt (t) — e At Tto the transition probability density of the Brownian particle 
velocity, conditioned by the initial data H*o at to — reads 



P(u '* |u ° )= ( ^(l- e -^ eXP ^ 2^(l-e-^) J ' (15) 

The process is Markovian and time-homogeneous, hence the above formula can be trivially ex- 
tended to encompass the case of arbitrary to 7^ : P (Tt, t\Tt o, to) arises by substituting everywhere 
t — to instead of t. 

Physical arguments (cf. demand (ii) preceding Eq. (4)) refer to an asymptotic probability distri- 
bution (invariant measure density) P(u) of the random variable Tt in the Maxwell-Boltzmann form 



p ™ ^{-^br) ■ (16 » 

This time-independent probability density together with the time-homogeneous transition density 
(15) uniquely determine a stationary Markovian stochastic process for which we can evaluate various 
mean values. Expectation values of velocity components vanish: (ui (t)) = UiP (Tt) dTt = for 
i = 1,2,3. The matrix of the second moments (velocity auto-correlation functions) reads 



/oo 
u i u j P(Tjt,t;Tto,to)dTtdTto (17) 
-oo 

where i,j — 1,2,3 and in view of P (Tt, t; Tto, to) = P (^t , t\Tt , to) P (Tt o) we arrive at the 
compact expression 



m 



-A\t-t \ 



m 



-P\t-t \ 



coso> c |t — to | sino> c |t — to | 
— sinu c |t — t 1 cosw c |t — t 
1 



(18) 
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In particular, the auto-correlation function (second moment) of the x-component of velocity equals 

(ui (t) mi (t )) = -^—e^ 1 *"' 01 C0SUJ c \t - to | (19) 

in agreement with white noise calculations of Refs. [Q and ||, cf. Chap. 11, formula (11.25). 
The so-called running diffusion coefficient arises here via straightforward integration of the function 
Rii{t) —< ui(t)ui(t ) > where r = t — t > 0: 

Dl (t) = f < Ul (Q)Mr) >dr = ^ + ^sin{w c t)- Pco S {u c t)]e X p(~pt) 
Jo m (V + u c A 

with an obvious asymptotics (the same for D 2 (t)): Dg = Zim t ^ oc Z? 1 (t) = ~~ 2 an d the large 
friction (lo c fixed and bounded) version D = ^jS-. 

3 Spatial process - dynamics in the plane 

The cylindrical symmetry of the problem allows us to consider separately processes running on the 
XY plane and along the Z-axis (where the free Brownian motion takes place). We shall confine 
further attention to the two-dimensional Xy-plane problem. Henceforth, each vector will carry two 
components which correspond to the x and y coordinates respectively. We will directly refer to the 
vector and matrix quantities introduced in the previous section, but while keeping the same notation, 
we shall simply disregard their z-coordinate contributions. 

We define the spatial displacement T* of the Brownian particle as folows 



— T^o = / ~tt (rf) dn (21) 
Jo 

where T? (t) is given by Eq. (^) (except for disregarding the third coordinate). 

Our aim is to derive the probability distribution of at time t provided that the particle position 
and velocity were equal T^o arL< i ^0 respectively, at time to = 0. To that end we shall mimic 
procedures of the previous section. In view of: 



r 



e-^lfo 



dse 



-Hv-s) 



A (a) 



(22) 



we have 
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is the inverse of the matrix A (regarded as a rank two sub-matrix of that originally introduced in 
Eq. (3)). Let us define two auxiliary matrices 



n = A-^l-e"**) (25) 
<j)(s) = A" 1 (l-e A ( s -*)) . 

Because of: 



e 



-At 



cxp^-t( P -% c )}=e-*( COSUJc \ SinLJ i )=e-^U(t) (26) 



uj c (3 J J V — sin uj c t cos u> c t 

we can represent matrices SI, (j) (s) in more detailed form. We have: 



_ _J [( P ^c\_- f 3t( P LUc \ f COSUJ ct sinu> c t i i 

^+^ c \\-uJ c J \ -w c /? ^ -sin^ c t cosw c t / ' 

and 



(s) = A' 1 (l - e-W-'^U (t - *)) = (28) 

_L / (3 lu c \ ( I -e^-V cos u c (s-t) -e^-VsmuJcis-t) \ 

0i + w 2 ^ /? J ^ e^( s -*) sinw c (s-t) 1 - e^( s -*) cos uj c (s - t) ) ' 

Next steps imitate procedures of the previous section. Thus, we seek for the probability distribution 
of the random (planar) vector it = f*<j) (s) (s) ds where it = — T^o — fiHV 

Dividing the time interval (0, t) into small subintervals to assure that (f> (s) can be regarded constant 
over the time span (jAt, (j + 1) At) and equal <f) (jAt), we obtain 

N-l r(j+l)At _ N-l N-l 

n = * (J At ) / ^ (s) d S = Y;<t> (ja*) ^ (At) = y, ^ ( 29 ) 

3=0 JjAt 3=0 j=0 

where T*,- = (jAt) £? (At) = ^ (At). 

By invoking the probability distribution (10) we perform an appropriate change of variables: ~r* j = 
4>j it (At) to yield a probability distribution of T*,- 

55 [T*j] - det «, [^r^] - ^-L-ttf . (30) 

Presently (not to be confused with previous steps (11)-(15)) we have 

det (s) = -g^-ji (l + e 2/3(s -* } - 2e*-'> cos lu c (s - t)) (31) 



G 



and 



r 1 w = Y 



+ e 2f)(s-t) _ 2e P(s-t) CQSUJc ( s _ t ) 

So, the inverse of the matrix <fij has the form: 



l_ e /»(»-t)tf(_( s _ t )) 



A . 



where 



and 



7j 



1 - e 0O' A *-*) cos lo c {jAt - t) e W At "*) sinw c (jAt - t) 



/3 -w c 



7j 



x + e 2/30A t -t) _ 2e /JWAt-t) CQSWc y At _ f) _ 



There holds: 



det ( 



7i 



and as a consequence we arrive at the following probability distribution of 



w[r j\ = 



1 / 1 



exp 



A, 



2 \ 



7] 4<7 At 



In view of 



Aj ( Jx 



(/3 2 +^) 7j (rj) 2 +(rjf) 



that finally leads to 



w r 



\AnqAtjjJ \ iqAtjj 



Since this probability distribution is Gaussian with mean zero and variance crj = 2qAt 
the random vector it as a sum of independent random variables T% has the distribution 



(*) 




*«P|-*Wl • (40) 



* 2 =£^=^£^ * 7 , ■ (41) 

i i c 

In the limit of At ->0we arrive at the integral 



ff2 = V"bii ^ s)ds (42) 

with J* 7 (s) ds = t + 0, where 

e = €) ( f ) = 4 (1 - e ' 2/jt ) - 2 o2 1 2 [/3+(^sin^-/?cosc^)e-' 3t ] . (43) 

That gives rise to an ultimate form of the transition probability density of the spatial displacement 
process: 

P (T, t |r„. - 0. *„) - - - - e.p (- Lg-g-^ j (44, 

with f2 = fi(t) defined in Eqs. (|25|), (27). Notice that an asymptotic diffusion coefficient Db = 
D pf + u* encodes an attenuation signature for the spatial dispersion (when uj c grows up at {3 fixed). 

The spatial displacement process governed by the above transition probability density surely is not 
Markovian. That can be checked by inspection: the Chapman-Kolmogorov identity is not valid, like 
in the standard free Brownian motion example where the Ornstein-Uhlenbeck process induced (sole) 
spatial dynamics is non-Markovian as well. 

4 Phase-space process 
4.1 Axial direction 

We take advantage of the cylindrical symmetry of our problem, and consider separately the (free) 
Brownian dynamics in the direction parallel to the magnetic field vector, e.g. along the Z-axis. 

That amounts to a familiar Ornstein-Uhlenbeck process in its extended phase-space form. In the 
absence of external forces, the kinetic (Kramers-Fokker-Planck equation) reads: 

dtW + uV z W = pVu(Wu)+qA u W (45) 

where q — D/3 2 . Here (3 is the friction coefficient, D will be identified later with the spatial diffusion 
constant, and (as before) we set D — ksT /m(3 in conformity with the Einstein fluctuation-dissipation 
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identity. The joint probability distribution (in fact, density) W(z,u, t) for a freely moving Brownian 
particle which at t = initiates its motion at xq with an arbitrary inital velocity uq can be given in 
the form of the maximally symmetric displacement probability law: 



W(z,u,t) = W(R, S) = [4n 2 (FG ■ exp{- ^^±^ } (46) 

where R = z - u (l - e^)/?" 1 , S = u- u e~^ while F = %(20t - 3 + 4e~^ - e~ 2 ^) G = 
D/3{1 ~ e" 2 ' 34 ) and H = D(l - e"' 3 *) 2 . 

4.2 Planar process 

Now we shall consider Brownian dynamics in the direction perpendicular to the magnetic field if, 
hence (while in terms of configuration-space variables) we address an issue of the planar dynamics. We 
are interested in the complete phase-space process, hence we need to specify the transition probability 
density P (~r* , if, t\~r*o, if o, io = 0) of the Markov process conditioned by the initial data if = if o 
and = T^o at time to — 0. That is equivalent to deducing the joint probability distribution 
of random vectors if and it, previously defined to appear in the form if = if (t)—e A4 lf o 
and it = — T^ — filf o respectively. Let us stress that presently, all vectors are regarded as two- 
dimensional versions (the third component being simply disregarded) of the original random variables 
we have discussed so far in Sections 2 and 3. 

Vectors if and R both share a Gaussian distribution with mean zero. Consequently, the joint 
distribution W {s, r) is determined by the matrix of variances and covariances: C — (cy) = ((xiXj)), 
where we abbreviate four phase-space variables in a single notion of x = (Si, S2, Ri, R2) and label 
components of x by i, j = 1, 2, 3, 4. In terms of it and if the covariance matrix C reads: 



/ (S1S1) (SiS 2 ) (S1R1) (S1R2) \ 
„ _ (S2S1) (S2S2) (S2R1) (S2R2) 
(R1S1) (RiS 2 ) (R1R1) {R1R2) 

V (r 2 Si) (R2S2) (r 2 Ri) (R2R2) J 

The joint probability distribution of if and R is given by 



W {tit) =W(1?) = ^ (_L_j " cxp ^Ec-^J (48) 
where c^Menotes the component of the inverse matrix C _1 . 

The probability distributions of if and it, which were established in the previous sections, deter- 
mine a number of expectation values: 

.9 = = (S2S2) = ^ (1 - e-^) (49) 

m 

while (S V S 2 ) = (S 2 Si) = 0. Furthermore: 



(47) 
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/ = (RiRj) = (R2R2) = 2 ksT M P „ (t + 6) = 2£> B (f + 6) . 



(50) 



m (3 2 + lu 2 

In addition we have (R1R2) = (R2R1) =0. As a consequence, we are left with only four non- 
vanishing components of the covariance matrix C: C13 = c 31 = (S1R1), c 14 = c 41 = (Sii? 2 }, c 23 = 
C32 = (S2R1), C24 = C42 = (S2R2) which need a closer examination. 

We can obtain those covariances by exploiting a dependence of the random quantities 5* and it 
on the white-noise term A (s) whose statistical properties are known. There follows: 



Si= [ dse~ 0{t ~ s) [cosui c (t-s)A 1 (s)+smuj c (t-s)A2(s)] 
Jo 



(51) 



S 2 = [ dse~ /3(t ~ s) [-sinu c (t- s) Ax (s) + cos uj c (t- s)A 2 (s)} 
Jo 



Ri 



I ds 7a\ 2 \P f 1 ~ e- fi{t - s) cosu c (t - 3)) + co.e-^-^ smu c {t - a)] A 1 {a) + 

fa-Pits) s - muJc _ s ) + bJc fi _ e -p(f-a) CQSU)c ^ _ s ^ A<1 ( s ) 



da— 5 5 

/3 2 + u 2 c 



Ri 



ds- 



1 



P 2 +" 2 

Io ds WT^ 



(l - e- f3{t - s) cosw c (t - s)) + f3e- l3{t - s) sinw c (t - *)] Ai (s) + 
(i-s)+/3(l-e-"<*- a >cosw c (t-*)) A 2 (s) . 



w c e ^* s) sinw c 



Multiplying together suitable components of vectors and it and taking averages of those prod- 
ucts in conformity with the rules {At (a)) = and {At (a) Aj (s 1 )} = 2qSij5 (s — s 1 ), where g = 
i,j = 1,2, 3, we arrive at: 



1 rt 

h={R 1 S 1 ) = {R 2 S2) = 2q— . / dse- f3( - t ~ s \(3cosuj c (t-s) + 

P + w c JO 



(52) 



and 



1 



w c sin lo c (t - a) - Pe-^-^} = q—„ 5 (l - 2e" /3t cos w c t + e" 2 ' 3 *) 



1 /•* 

k = (iixSa) = - (i? 2 5i) = 2g-= . / dae-^-^I-^sinwe (t - a) + 

P + w c JO 



(53) 



iv c cos w c (i - s) - iv c e t3(t s) ] = q—^- — ^ 



2e-f }t smu; c t-^(l~e- 2 e t ) 
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The covariance matrix C — (cjj) has thus the form 



h 



C 



while its inverse C reads as follows: 



9 



( 9 


h k f 
\ -k h 



-k \ 
h 


/ J 



(54) 



1 



detC 







o / 

-h -k 
\ k -h 



-h k \ 

-k -h 

9 

g J 



(55) 



where det C = (fg — h 2 — k 2 ) 2 . The joint probability distribution of "3? and it can be ultimately 
written in the form: 



W 



4tt 2 (fg -h 2 -k 2 ) 



exp 



5* 



R 



2h 5" • 7? + 2k ( 5" x 7?) ^ 



2 (/«? - ^ _ fc2) 



(56) 



(57) 



In the above, all vector entries are two-dimensional. The specific i = 3 vector product coordinate 
in the exponent is simply an abbreviation for the (ordinary i? 3 -vector product) procedure that involves 
merely first two components of three-dimensional vectors (the third is then arbitrary and irrelevant), 
hence effectively involves our two-dimensional it and 7?. 
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